Detection of a virial shock around the Coma galaxy cluster 
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Galaxy clusters, the largest gravitationally bound objects in the Universe, are thought to grow by 
accreting mass from their surroundings through large-scale virial shocks. Due to electron acceleration 
in such a shock, it should appear as a 7-ray, hard X-ray, and radio ring, elongated towards the large- 
scale filaments feeding the cluster, coincident with a cutoff in the thermal Sunyaev-Zel'dovich (SZ) 
signal. However, no such signature was found so far, and the very existence of cluster virial shocks 
has remained a theory. We report the discovery of a large, ~ 5 Mpc diameter 7-ray ring around 
the Coma cluster, elongated towards the large scale filament connecting Coma and Abell 1367. The 
7-ray ring correlates both with a synchrotron signal and with the SZ cutoff. The 7-ray, hard-X-ray, 
and radio signatures agree with analytic and numerical predictions, if the shock deposits a few 
percent of the thermal energy in relativistic electrons over a Hubble time, and ~ 1% of the energy 
in magnetic fields. The implied inverse-Compton and synchrotron cumulative emission from similar 
shocks dominates the diffuse extragalactic 7-ray and low frequency radio backgrounds. Our results 
reveal the prolate structure of the hot gas in Coma, the feeding pattern of the cluster, and properties 
of the surrounding large scale voids and filaments. The anticipated detection of such shocks around 
other clusters would provide a powerful new cosmological probe. 



1. INTRODUCTION AND OVERVIEW 

In the hierarchical paradigm of large scale structure 
(LSS) formation, galaxy clusters are the largest object 
ever to virializc. With masses in excess of 10 14 Mq, they 
are located at the nodes of the cosmic web, where they 
accrete matter from the surrounding voids and through 
large scale filaments. Due to their vast size, galaxy clus- 
ters resemble island Universes seen at great distance, pro- 
viding a powerful cosmological probe and a unique astro- 
physical laboratory. 

The gas accreted by a cluster is thought to abruptly 
heat and slow down in a strong virial shock wave 
surrounding the cluster. Such collisionless shocks are 
thought, by analogy with supernova remnant shocks, 
to accelerate charged particles to > TcV energies, 
where they Compton-scatter cosmic-microwave back- 
ground (CMB) photons up to the 7-ray band 1- ' 5 . Conse- 
quently, one expects to find 7-ray rings around clusters , 
as indicated by cosmological simulations ' , which sug- 
gest an elliptic morphology elongated towards the large- 
scale filaments feeding the cluster''. Such rings are also 
expected in hard X-rays' 1 , and should coincide with a syn- 
chrotron radio ring > and with a cutoff in the thermal 
Sunyaev-Zel'dovich (SZ) signal 8 . 

However, no such virial shock signature was detected so 
far, although a stacking analysis of EGRET data around 
a sample of 447 rich clusters did show a 3a signal 9 . Upper 
limit were imposed on the 7-ray emission from clusters 
such as Coma 10-13 , but mostly focusing on the central 
parts of the cluster, well within the virial radius. The 
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very existence of cluster-scale virial shocks has thus re- 
mained unconfirmed. 

The Coma cluster (Abell 1656) is one of the rich- 
est nearby clusters. With mass M ~ 1O 15 M0, tem- 
perature T ~ 8 keV, and richness class 2, it lies only 
~ 100 Mpc away 14 , at a redshift z = 0.023. Its virial 
radius, R ~ -R200 — 2.3 Mpc, corresponds to an angular 
radius vb ~ ^> 2 oo — 1-3°, where subscript 200 refers to 
an enclosed density 200 times above the critical value. 
The cluster is somewhat elongated in the East- West di- 
rection, in coincidence with the Western LSS filament 15 
that connects it with the cluster Abell 1367 (see Figure 
6a). There is X-ray, optical, weak lensing, radio 1 ' 1 , and 
SZ evidence that the cluster is accreting clumpy matter 
and experiencing weak shocks towards the filament well 
within the virial radius, at ij) ~ 0.5° radii. 

The VERITAS Cerenkov array has produced a d ~ 
5° diameter 7-ray mosaic 13 of Coma, at energies E > 
220 GcV. We argue that the significance map (Figure la) 
shows extended 7-ray emission away from the center, that 
appears as a thick (~ 0.5°) elliptical ring with semi-minor 
axis b ~ 1.3°, elongated along the East- West direction, 
with semi-major to semi-minor axes ratio £ = a/b > 3. 
The ring is seen at a 4.5a confidence level (effectively 
a single-parameter fit, estimated using mock VERITAS 
noise maps; see Figures 2 and 3 and §2). However, the 
data was taken in w = 0.5° wobble mode 18 , inappropriate 
for extended emission on larger scales; correcting for the 
excessive background removal indicates a > 5.7cr signal. 

X-ray and galaxy maps suggest that the LSS filament 
and the semi-major axis of the thermal gas distribution 
are slightly tilted Southwest, at an angle — 10° < <j> S 0°. 
Indeed, the significance of the 7-ray structure is maxi- 
mized when thus tilted; the best fit is obtained for ( ~ 4.5 
and (f> ~ —5°. If no assumption is made regarding 0, the 
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(a) VERITAS significance map (b) Simulated cluster in 7-rays 

FIG. 1: Observed and simulated 7-ray maps of Coma. 

Left: VERITAS > 220 GeV significance map 1:i of Coma for 9 = 0.2° integration (illustrated by the red cir- 
cle). Elliptical bins are shown (thin dashed cyan contours) for Ab = 0.2°, £ = a/b = 2 and <fr = —5°. The 
bins showing enhanced emission are highlighted (bounded by thick, long-dashed green curves). Right: Simulated 
map of a Coma-like cluster from a ACDM simulation 5 . The largest in the 200 Mpc simulation box, this cluster 
has' mass M ~ 10 15 Af Q and temperature T ~ 8 kcV. The 8.5° diameter image was convolved with a ~ 0.23° 
beam, comparable to the VERITAS map, and rotated such that the large-scale filament extends to the West. 
Colorbar: log 10 (J/10 -8 cm -2 s _1 ster -1 ) brightness above 220 GeV for u£ e = 5%. The regions corresponding to 
the VERITAS mosaic (solid cyan contour) and to the VERITAS ring (elliptic dashed contours) are highlighted. 



ring significance weakens to 4.2(7 (> 5.4(7 corrected). 

The size and morphology of the VERITAS ring re- 
semble the predicted 7-ray signature of a Coma-like, 
M ~ 1O 15 M and T ~ 8 keV cluster in a ACDM 
simulation'. The simulated map (Figure lb), oriented 
such that the major simulated LSS filament lies to the 
West and normalized for a nominal 3.3a wobble-mode 
VERITAS detection (no remaining free parameters) , cor- 
relates with the VERITAS map at the 3.7a (4.9a cor- 
rected) confidence level. 

The ring brightness agrees with predictions too, if the 
relativistic electron fraction £ e of the thermal energy ac- 
cretion rate u (normalized by the total thermal energy 
over a Hubble time) is a few percent. The associated 
~ 5 eV s _1 cm~ 2 flux at 20-80 keV energies can ex- 
plain the hard X-ray signal observed in Coma 19 . Sim- 
ilar, £ e M ~ 4% values are inferred from the 3c EGRET 
stacking analysis' ', and from LSS hard X-ray sources 6,20 . 

In order to examine the coincident radio signals, we 
spectrally decompose the WMAP seven-year, full-sky 
maps" into synchrotron and thermal SZ maps. Exclud- 
ing the inner, b < 1° part of the cluster, the VERITAS 
data correlates with the synchrotron map at a +2.8(7 
(+3.9cr) confidence level (Figure 9), and anti-correlates 



with the SZ signal at — 2.6a (— 3.6cr). The radio signals 
are stronger in the Western half of the ring, towards the 
LSS filament, where the correlations strengthen up to 
+3.5ct (+4.6<7) with the synchrotron and — 3.9a (—5.0a) 
with the SZ. 

These peripheral radio signals agree with predictions 
for the synchrotron emission from the 7-ray emitting 
electrons, and for the SZ decline with thermal pressure 
beyond the shock, provided that B ~ 0.6 fiG magnetic 
fields are found downstream, corresponding to £b — 1% 
magnetization. Bright spots along the ring suggest local 
enhancements, up to u£, e £,B ~ a few times its mean value. 

While the SZ signal, extending in some parts out to 
ip ~ 2.3°, agrees with SZ detections by WMAP 21 and 
Planck 1 ', the peripheral synchrotron signal coincident 
with the 7-rays and with the SZ cutoff was not reported 
so far. Note that the 7-ray ring and these radio features 
are much farther out than the tp ~ 0.5° weak shocks pre- 
viously reported based on X-ray, radio 11 ', and SZ 1 ' data. 

Inside the b = 1° ellipse, the VERITAS and syn- 
chrotron maps show a —1.9a (—2.7a) anti-correlation, 
mainly (—2.2a (—2.8a)) in the West, evident from the 7- 
ray undcrluminous radio halo and relic (Figure 9). Such a 
signature, predicted in secondary models where cosmic- 
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FIG. 2: Significance of the VERITAS ring. 

Elliptical binning (symbols; connecting lines are guides 
to the eye) of the significance map (Figure la) with 
A6 = 0.6°, an East- West orientation (<fi = 0), and 
semi- major to semi- minor axes ratios £ = a/b = 1 
(circles), 2 (diamonds), 3 (squares), 5 (right triangles) 
and 1000 (stars). A slightly tilted, (j> = —5° orienta- 
tion provides a slightly better fit, illustrated for ( = 5 
(down triangles and solid line). The yellow band cor- 
responds to the bright extended ring marked in Figure 
la. The detection significance of the ring can be es- 
timated from the enclosed flux (the ordinate gives the 
nominal significance S), but more accurate estimates are 
found by considering a large sample of mock maps con- 
taining only noise, and applying the same elliptical tem- 
plates, either constrained to have an East- West elonga- 
tion (labelled dashed lines) or picking the angle <fi that 
maximizes the significance in each map (labelled dot- 
dashed lines). In either case, a correction due to exces- 
sive wobble-mode background subtraction is necessary. 



ray ion collisions inject relativistic electrons in both 
strongly (> 3 [iG) and weakly magnetized regions, which 
consequently appear radio bright/7-ray faint and vice 
versa, is unnatural in primary models where electrons are 
(re) accelerated in turbulence 25 ' 21 or weak shocks. Sec- 
ondary electrons can explain radio halos ' , relics and 
minihalos ' as arising from the same cosmic-ray distribu- 
tion, if the latter is nearly homogeneous 22 . The extent 
and magnitude of the anti-correlation observed, if con- 
firmed, support such secondary models, and favor an ex- 
tended cosmic-ray distribution. 

A future VERITAS observation deeper than the 
present 18.6 hour exposure would decisively test our re- 
sults, in particular if wobble-mode corrections are mini- 
mized. This would allow a more precise reconstruction of 
the virial shock and a measure of variations in the elec- 
tron deposition rate along the shock front. Such data is 
valuable in the study of LSS formation at low redshifts, 
mapping the feeding pattern of Coma, probing the sur- 
rounding voids and filaments, and tracing the warm-hot 
intergalactic medium (WHIM) immediately behind the 



shock. 

A detailed, more significant anti-correlation between 
7-ray and synchrotron maps at inner, ip < 1° radii would 
confirm the hadronic, secondary electron model for clus- 
ter radio sources, and provide a direct measure of mag- 
netic variations in the ICM. A deeper 7-ray observation 
would thus probe the cosmic-ray proton distribution out 
to large radii, address the origin of these protons, and 
place an upper limit on proton acceleration in the virial 
shock. 

The VERITAS energy threshold, E ~ 220 GeV, is 
not far below the expected photon cutoff, E max ~ 
50(S/0.6 fiG) TeV. Observations at higher energies 
would become increasingly sensitive to this cutoff, thus 
probing the conditions at the shock. Interestingly, no ev- 
idence for extended emission was found in a 1.1 TeV ob- 
servation of Coma by the HESS telescopes 28 , with ~ 40% 
the sensitivity in the VERITAS mosaic, but the 7-ray 
ring may be shallowly buried in the HESS noise. 

Calibrating analytic 1 '', and numerical''' estimates 
with the above w£ e and £b values, yields E 2 dJ/dE ~ 
0.3 keV(£ e u/5%) s _1 cm -2 stcr -1 inverse- Compton and 
vl v ~ 3 x 10- 11 (^ e u/5%)(Cs/l%) erg s" 1 cm" 2 stcr" 1 
synchrotron cumulative backgrounds from cluster shocks. 
These constitute fair fractions of the 7-ray and low- 
frequency radio signals observed far from the Galactic 
plane, and dominate the respective diffuse extragalactic 
backgrounds ' . 

The paper is organized as follows. In §2 we show 
that an elongated 7-ray ring around Coma was observed 
by the VERITAS Cercnkov array. Methods to correct 
for excessive wobble-mode background subtraction are 
discussed in general in §3, and specifically for cross- 
correlation analyses in §4. Contamination by Galactic 
foregrounds is shown to be negligible in §5. The 7-ray 
ring is shown to agree with simulation-based predictions 
in §6. Significant correlations with both synchrotron and 
SZ maps are presented in §7. In §8 we show that the 7- 
ray signal agrees with predictions, if the shock deposits a 
few percent of the thermal energy per Hubble time in rel- 
ativistic electrons. In §9 we show that the synchrotron 
and SZ signals also agree with the predictions, if mag- 
netic fields carry £g — 1% of the thermal energy. A 
7-ray-radio anti-correlation inward of the 7-ray ring is 
interpreted in §10 as possible preliminary evidence for 
strong local magnetization and secondary electron injec- 
tion at small radii. The results are summarized and dis- 
cussed in §11. A simple j3- model of the shock is provided 
in Appendix A. 

We adopt a concordance flat ACDM model with nub- 
ble constant H = 70 km Mpc" 1 , a baryon fraction fi, = 
17%, and a hydrogen mass fraction \ = 0.75. The plasma 
is approximated as an ideal gas with adiabatic index 
r = 5/3 and mean particle mass /i = 0.59m p , where 
m p is the proton mass. Colorbars are square- weighted. 
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2. VERITAS 7- RAY RING AROUND COMA 

VERITAS has produced a d ~ 5° diameter 7-ray 
mosaic 1 ' 5 of Coma, at energies > 220 GeV. The signifi- 
cance map (Figure la), taken in w = 0.5° wobble mode 
with a 9 = 0.2° integration radius, shows some extended 
7-ray emission away from the center. The 7-ray struc- 
ture appears as an elliptical ring, elongated along the 
East- West direction, or as two parallel filaments lying 
symmetrically both North and South of the cluster. The 
spectrum in the VERITAS band is probably flat, as a 
p = 2.4 photon spectral index was used to optimize the 
gamma-hadron separation cuts . In §8 we show that a 
flat spectrum with index p = 2 (equal energy per loga- 
rithmic energy interval) is consistent with observations 
at lower energies. 

Denote the VERITAS significance map by Sj. The 
data is taken in wobble mode 18 , where a pixel of signifi- 
cance Sj means that in a beam of radius 9 = 0.2° centered 
upon the pixel, there are SjN a counts in excess of the an- 

1 12 

ticipated A0.2 background events, where N a = N 2 and 
the background is estimated based on a w ~ 0.5° ring 
around the pixel. A 8 2 analysis 1 ' shows a background of 
about iVo.2 — 1000 photon counts above 220 GeV in such 
a = 0.2° bin. The background per pixel is denoted by 
Ni, = Nq. 2/^0.2, where no. 2 is the number of pixels in the 
beam. First we assume that the background is uniform, 
so the pixels are not correlated before the beam integra- 
tion. The nominal significance of an extended source may 
then be estimated as 

s= aN ^=-^—, (1) 

where s = s j an d n arc respectively the total signifi- 
cance and the number of pixels enclosed. 

In order to quantify the 7-ray structure and assess its 
significance, we fit the data with a thick, elliptical ring 
model. The center of the ellipse is chosen as the ROSAT 
X-ray peak. The median semi-major axis a is taken along 
the East-West direction, as inferred from the ROSAT 
map and from the orientation of the SDSS galaxy fila- 
ment. The thickness of the ring is fixed at Ab = 0.6° 
along the semi-minor axis, because a thicker ring would 
exceed the 2w = 1° limit imposed by the wobble-mode 
background subtraction 5 ". The ring is then defined by its 
median semi- minor axis b and elongation ratio £ = a/b, 
taken as two free parameters. However, for sufficient 
elongation, £ > 3, the results (Figures 2 and 3) de- 
pend weakly upon £. The (effectively single-parameter) 
fit then indicates an extended structure around b ~ 1.3°, 
found at an S = 3.3er nominal significance level. 

The SDSS (Figure 6a) and ROSAT maps suggest that 
the LSS filament and the semi-major axis of the gas dis- 
tribution are slightly tilted towards the Southwest, at an 
angle —10° < <j> < 0°. Indeed, the significance of the 7- 
ray structure is maximized when tilted at precisely such 
angles (Figure 3). The best fit is obtained for ( ~ 4.5 
and </> ~ -5°. 
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FIG. 3: Orientation of the VERITAS ring. 

The nominal significance S (thin, solid curves) and cor- 
responding semi-minor axes b (thick, dashed curves) 
are shown for elongated rings of various ellipticity 
C (legend). The yellow sector corresponds to the 
-10° < 4> < 0° orientation of the LSS filament. 
The secondary peak towards <fi ~ 70° with b ~ 
0.5° is related to the primary peak: the semi-major 
axis of the small ellipse solution roughly coincides 
with the semi-minor axis of the large ellipse solution. 

A more reliable way to estimate the significance of the 
extended emission is to apply the same elliptical tem- 
plates to a large sample (~ 1000) of mock significance 
maps, containing random noise but no signal. We prepare 
these maps by injecting Poisson noise at the observed 
background level to the VERITAS mosaic template, in- 
tegrating over the same 9 = 0.2° beam, and converting 
to a significance map using either the mean background 
or a simulated w = 0.5° wobble-mode background. 

None of the resulting mock maps shows a signal as 
strong as the observed S = 3.3<r, so we extrapolate to 
the high significance levels assuming a Gaussian distri- 
bution. The resulting significance levels are shown as 
dashed curves in Figure 2, for elliptical rings elongated in 
the East- West direction (based on 1000 maps at a given 
£) . If we make no a-priori assumption regarding the ring 
orientation 0, but rather choose in each mock map the 
value <j) that maximizes its nominal ring significance, we 
obtain higher mock detection levels, shown as dot-dashed 
contours in the figure (based on 200 maps, each scanned 
with Acf> = 20° tilt intervals) . The results depend weakly 
(within the statistical fluctuations) on £. 

The mock statistics indicate that the VERITAS ring 
presents at the 4.5er confidence level if we assume that it 
is aligned with the LSS filament, and at 4.2a confidence if 
we make no a-priori assumption regarding <j). These con- 
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servative estimates are based on the mean background, 
and are approximately independent of £; higher, ~ 5.5cr 
significance levels are obtained when wobble-mode is used 
to estimate the mock background. As we show in §3 be- 
low, the detection significance for such extended struc- 
tures is actually higher, > 5.7a for East- West alignment, 
and > 5.4 for arbitrary </>. 

The signal is locally on the order of < 10% of the 
noise, and stronger than the average Galactic foreground 
at these energies 31 (see §8). In the Coma region, the 
expected Galactic foreground is relatively low, spatially 
smooth, and featureless on < 0.5° scales; see §5. Wobble- 
background subtraction is highly efficient in removing 
smooth foregrounds; a disadvantage is that it also re- 
moves much of the signal for extended sources. 

In wobble mode, an estimated background is sub- 
tracted based on a ring of radius w around each pixel. 
This method, appropriate for point sources, distorts and 
smears out features extended over scales comparable or 
larger than w. As we show below, analyses of the 7-ray 
features, of a histogram of the 7-ray pixels, and of wobble 
simulations, indicate that the wobble subtraction reduces 
the significance of the extended features observed by at 
least AS = 1.2a. Therefore, the corrected significance of 
a 7-ray ring aligned with the LSS filament is estimated 
as > 5.7a. 



3. CORRECTING FOR WOBBLE-MODE 
BACKGROUND SUBTRACTION 

Consider the wobble-induced correlation between 7-ray 
pixels. Due to the wobble background subtraction, the 
significance of a pixel is anti-correlated with the pixels 
w ~ 0.5° around it. One implication is that sources 
that extend over scales much larger than w are effec- 
tively erased by the background subtraction, except in 
an ~ w wide band in their periphery. In such a band, 
and in sources that extend < w in at least one direction, 
the signal is diminished but not completely erased. For 
example, the estimated signal at the inner edge of a large 
linear source will be diminished by a factor of two, and 
the signal at the outer edge will have an equal but nega- 
tive significance. Therefore, for sources that extend over 
> w in at least one direction, the actual significance may 
be > 2 times higher than estimated from the wobble- 
mode map. 

In order to quantify the effect wobble-mode has on 
the signal, we examine the distribution of 7-ray signifi- 
cance pixels, using the same 9 = 0.2° integration beam 
applied to the map. The significance histogram (Figure 
4), taken with 5a = 0.1 intervals, can be modeled 13 as 
a single Gaussian function, with adjusted coefficient of 
determination R 2 = 0.9956. 

However, the histogram shows peak deviations and 
sharp cutoffs, indicating a more complicated underly- 
ing distribution. Indeed, a simple model combining two 
Gaussian functions provides a much better fit to the data 



(R 2 = 0.9967; note that R improves slowly near 1). The 
best fit is obtained for Gaussian means fii = — 1.23±0.21 
and H2 = 0.23±0.15, with corresponding standard devia- 
tions 0.69±0.12 and 0.91±0.08. The data were also mod- 
elled as the sum of two Gaussian functions constrained to 
have the same standard deviation; the resulting best fit 
is similar to the unconstrained sum: \i\ = —0.96 ± 0.06 
and jU2 = 0.41 ±0.05, with standard deviation 0.83 ±0.04 
(R 2 = 0.9967). 

Both fits suggest that a signal is present in the map, 
with pixels that are on average [12 — \i\ — lAa above the 
noise. Interestingly, a histogram of the pixels found only 
in the central, ~ 2° wide region of the map agrees much 
better with a single Gaussian !_> . This indicates that the 
additional signal is found at radii > 1°, in agreement 
with the inferred ring. 

The above arguments indicate that the local pixel sig- 
nificance should in fact be Asj ~ — fXi ~ 1.2<7 higher 
than it appears in the significance map (Figure la). This 
imposes a lower limit AS = 1.2(n/no.2) 1 ^ 2 o' to the cor- 
rection factor needed in order to estimate the significance 
of an extended structure or of a cross-correlation signal 
(see §4) with respect to the true background. We conser- 
vatively use AS* = 1.2cr for extended emission exceeding 
the beam size. 

Another indication that such a correction is needed 
are the negative, < — 3cr significance regions, lying just 
outside the bright feature we identify as a ring. Assum- 
ing that the bright region is part of an extended, linear 
structure, this suggests that its local significance should 
be Asj > (3 + £ii)/2 ~ 0.9<7 higher than it appears. 

Another way to assess the effect of wobble background 
subtraction is to apply it to a simulated cluster map. We 
thus normalize the 7-ray map of the simulated cluster de- 
rived in §6 below, add homogeneous poisson noise, intro- 
duce a w = 0.5° wobble-mode correction, and integrate 
over the 9 = 0.2° beam. The normalization and back- 
ground are chosen such that the final map thus obtained 
shows 4.0cr bright spots over a background comparable to 
the VERITAS background. For such parameters, a map 
identically prepared but with mean background, rather 
than wobble background subtraction, shows 6.1cr bright 
spots, implying a true significance ~ 50% higher than 
inferred in wobble-mode. 

A 7-ray structure that extends over scales much greater 
than w would be largely erased by the wobble background 
subtraction. Hence, the above estimates provide only 
a lower limit to the correction needed to determine the 
significance of such a structure. 



4. CROSS-CORRELATION SIGNALS AND 
WOBBLE-MODE CALIBRATION 

Next, consider a cross-correlation between the VERI- 
TAS map Sj and some potential tracer of the virial shock, 
qj . After binning the two maps onto some common grid, 
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FIG. 4: Histogram of pixel significance in the VERITAS map. 

A single Gaussian function (dot-dashed) provides a marginal fit to the data. A much better fit is provided by the 
sum (solid) of two Gaussian functions (dashed); see text for function parameters. Inset: logarithmic histogram. 



we compute the correlation as 

n _ ST (sj - Ms) (gj - Mg) 



where the means fi, standard deviations a, and number 
of pixels n, pertain to the region being examined. The 
factor no. 2 accounts for the correlation between pixels 
within the 8 = 0.2° beam, ensuring that the result is 
approximately independent of the grid resolution. The 
significance of the correlation is determined by generat- 
ing a large (> n 2 , converged) sample of maps with ran- 
domly swapped pixels, and computing the corresponding 
significance distribution. 

Consider the effect wobble-mode background subtrac- 
tion has on such cross-correlations. As in the significance 
estimation discussed in §3, any existing correlation is di- 
minished because extended bright regions experience ex- 
cessive background subtraction. However, the effect here 
is more severe, because background regions lying < w 
away from the bright structure are affected as well, ap- 
pearing excessively faint. As such "dark edges" appear 
in the wobble-mode map but not in the tracers, the cor- 
relation is further diminished. 

The effect of wobble-mode background subtraction on 
a cross-correlations signal can be crudely estimated by 
correlating mock tracers with simulated 7-ray maps, pre- 
pared with (e.g., Figure 7) and without (e.g., Figure lb) 
wobble background subtraction. We generate mock trac- 
ers by adding random noise of various mean amplitude to 
the significance map of the simulated cluster, smoothed 
with a Gaussian filter of standard deviation r q . The re- 
sults are well-fit by a wobble correction factor 

AC~0.18e (r * /a2O) (d/5°)C , (3) 

where C is the nominal significance of the cross- 
correlation using the wobble-mode signal, C + AC is the 



true significance, and 7r(d/2) 2 is the solid angle involved. 

For example, for the typical r q = 0.25° of WMAP, the 
C ~ +2.8cr correlation found (see §7 below) between 7- 
ray and synchrotron emission outside the inner 6=1° 
ellipse bounding the 7-ray structure from below yields a 
corrected C + AC = +3.9er confidence level. As pointed 
out in §3, these corrections arc only lower limits if the 
7-ray structure is more extended than the simulated ring 
used to derive Eq. (3). 

5. NEGLIGIBLE GALACTIC CONTAMINATION 

When averaged over large angular scales, the Galactic 
foreground anticipated above 220 GeV based on Fermi 
observations ' 1 is lower by at least a factor of a few than 
the ring signal (see §8). Moreover, Galactic signals on 
> w = 0.5° scales are efficiently removed from the VER- 
ITAS map by the wobble-mode background correction. 
Next, we show that no significant Galactic contamina- 
tion on < w scales is expected in the VERITAS mosaic 
of Coma. 

The main 7-ray Galactic foregrounds at interme- 
diate latitudes arise from (i) inverse Compton emis- 
sion from cosmic-ray (CR) electrons; (ii) relativistic 
bremsstrahlung of CR electrons with the interstellar 
plasma; and (Hi) nucleon-nucleon scattering between CR 
ions and the interstellar medium. The CRs are unlikely 
to show significant structure on < 0.5° scales, as this 
would require < 10 pc magnetic fields at < 1 kpc dis- 
tances. Moreover, synchrotron maps ' show relatively 
very little emission from the Coma region. Therefore, 
the Galactic foreground on small scales is dominated by 
electron bremsstrahlung and by ir° decay from nucleoli 
collisions, both approximately following the distribution 
of gas along the line of sight. 

Hence, we may estimate the Galactic contamination 
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FIG. 5: Tracers of Galactic foreground near Coma. 

The Ha full-sky map ; ' (left; colorbar units: Raylcigh) and the SFD dust map ;: ' (right; color- 
bar: E(V-B) reddening) show only little Galactic structure in the region of the VERITAS mo- 
saic (enclosed by cyan contour), and in particular towards the 7-ray ring (dashed green contours). 



remaining in the VERITAS map using gas tracers, such 
as the full-Sky Ha 55 and infrared dust emission" 1 maps 
(Figure 5). Both tracers indicate a very low level of 
Galactic structure in the Coma region, and towards the 
7-ray ring in particular, with only minor structure found 
mainly in the East. We confirm this using a spectral lin- 
ear regression of the seven-year WMAP data, outlined 
in §7. Both tracers show no positive correlation with 
the VERITAS map, ruling out any significant remaining 
Galactic contamination. 

Note that synchrotron tracers should not be used here 
to estimate or eliminate the Galactic foreground. First, 
these tracers are inefficient as they follow the smooth CR 
electron distribution. More importantly, they include an 
interesting synchrotron signal from the Coma virial ring, 
as we argue in §7. 

6. GAMMA-RAYS FROM A SIMULATED 
CLUSTER 

The VERITAS signal can be assessed by comparing it 
to the 7-ray signature of a simulated Coma-like cluster in 
a ACDM simulation ! . The 7-ray signal was computed by 
injecting relativistic electrons at the strong shocks of the 
simulation, with a Mach number-dependent power-law 
spectrum given by the Fermi diffusive shock acceleration 
model. These electrons were assumed to carry a fraction 
£ e of the thermal post-shock energy; see §8 for details. We 
choose the richest cluster in the (200 Mpc) 3 simulation 



box, found 7 to have a Coma-like mass ~ 1O 15 M and 
temperature T ~ 8 keV. 

The simulated virial shock (Figure lb) has an elliptic 
structure, elongated towards the main LSS filament (Fig- 
ure 6a), with semi-minor axis b ~ 1.3° and semi- major 
axis a ~ 2.5° (at the distance of Coma); the average 
< 220 GeV spectral index is p ~ 2.05. The shock elonga- 
tion, the morphology of the surrounding filaments, and 
the enhanced accretion towards the main filament, sug- 
gest a prolate geometry aligned with the main filament 
feeding the cluster. Unrelaxed, merger clusters such as 
Coma are indeed known to preferentially be prolate 5 ' . 

The center of the simulated cluster is determined by 
the peak density. The simulated map is rotated about 
this point such that the main simulated LSS filament 
coincides with the </> ~ —5° orientation of the SDSS fila- 
ment and the ROSAT ellipse. By coincidence, the agree- 
ment between the simulated filaments and the SDSS map 
remains fairly good even far from the cluster, if the sim- 
ulated map is inverted along the semi-major axis of the 
virial shock (Figure 6). 

The simulated 7-ray signal shows a positive correla- 
tion with the VERITAS map (Figure 1), at the 3.7a 
confidence level (4.9<r after correcting for the excessive 
wobble- mode background subtraction; see §3). The cor- 
relation remains strong, > 3.5ct (> 4.7er), if the maps 
are rotated by < 5° with respect to each other. The ex- 
tent of the virial shock may be somewhat overestimated 
by the adiabatic, smooth particle hydrodynamics (SPH) 
simulation. Indeed, the correlation strengthens to 4.1er 
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(a) SDSS galaxies with simulated LSS filaments (b) Corresponding simulated 7-ray map 

FIG. 6: The LSS environment of Coma. 

Left: SDSS galaxies at redshift 0.018 < z < 0.028 near Coma. The simulation roughly reproduces the filamentary 
structure near Coma, especially when reflected about the semi-major axis of the virial shock (solid red contours show 
the density- weighted T = 1 kcV gas in a 4.4 Mpc thick slice through the simulation box 3 ). Right: 7-ray map for 
£ e m = 5%, showing the > 220 GeV brightness as log 10 ( J/10~ 8 cm~ 2 s _1 stcr -1 ) after said reflection. Elliptic dashed 
contours from Figure la are superimposed, and the region shown in the right panel is marked by a dot-dashed circle. 



(5.3<r) if the simulated map is resized by a factor 0.9. 

The mass accretion rate M of the cluster, normalized 
by the mass-to- Hubble time ratio in the dimensionless pa- 
rameter m = M/(MH), varies among clusters and fluc- 
tuates in time, and is found' to be somewhat low when 
averaged over the simulated cluster, m ~ 0.3. The 7-ray 
map shown in Figure lb corresponds to normalized aver- 
age accretion rate m = 1, and an acceleration efficiency 
£ e = 5%. Equivalently, we may examine the accretion 
rate of the total thermal energy Uth, using the dimen- 
sionless parameter u = Uth/{UthH). In the simulation, 
we find -5 '' ii ~ m, although slightly larger, ii = (5/3)m 
values arise in a simple isothermal sphere model 4 . 



7. SIGNIFICANT CROSS-CORRELATIONS 
BETWEEN 7-RAY AND RADIO SIGNALS 

In order to quantify the coincidence between 7-ray and 
radio signals, we analyze the WMAP seven-year, full- 
sky maps 21 . We use a simple spectral decomposition of 
the five WMAP bands into synchrotron and SZ maps, 
avoiding any spatial masking which may complicate the 
comparison to the VERITAS map. 

The synchrotron spectrum is approximated as a pure 
power-law, with brightness l v oc v a and spectral index 
a = —1.2, typical of cluster radio sources; the results are 
not sensitive to the precise spectrum for — 2 < a < — 1. 



The SZ signal is clear in the central parts of Coma, but 
is difficult to spectrally separate from the CMB fluctua- 
tions at large radii 21 . For simplicity, we leave these CMB 
fluctuations in the SZ map, thus somewhat diluting any 
correlations with the VERITAS map. 

In general, the WMAP data includes non-negligible 
foregrounds from Galactic synchrotron, dust and free- 
free emission , which are strong towards the Galactic 
plane and in principal should be removed. However, as 
discussed in §5, these signals are particularly weak in the 
Coma region (Figure 5), are spatially too smooth (in par- 
ticular the synchrotron foreground) to have a 7-ray coun- 
terpart surviving the wobble-mode background removal, 
and their tracers do not correlate with the VERITAS 
map. 

Moreover, the extragalactic synchrotron signal cannot 
be separated from the Galactic a ~ — 1 synchrotron, 
a ~ —0.7 dust, and a ~ —0.35 free-free signals, using 
spectral regression alone, in particular considering the 
varying synchrotron spectral index in the cluster radio 
sources. Therefore, we leave these faint and probably 
smooth Galactic contaminations in our synchrotron ex- 
tragalactic tracer. This only dilutes the correlations be- 
tween this tracer and the VERITAS map, as the Galactic 
signals show no positive correlations with the latter. 
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FIG. 7: Wobble-mode background- 
subtracted 7-ray map of the simulated cluster. 

The map corresponds to the VERITAS parameters: 
At = 18.6 hour integration, A = 4 x 10 8 cm 2 ef- 
fective area, and 90% Poisson noise. While the 
bright spots are ~ 3ct for £ e rh = 5%, here we use 
£ e m = 50% to also show the dimmer Western edge. 



Therefore, we fit the five WMAP channels by 

AT = T — T cmb (4) 



y 



v coth 



-4 



T syn 



-3.2 



where v = hv /kBT cm b is the dimensionless frequency and 
y is the Comptonization parameter. The low-frequency 
SZ temperature shift AT SZ = — 2y and the synchrotron 
brightness temperature T syn (at the arbitrary frequency 
vq) are used as fit parameters, giving the tracer maps 
(Figure 8). 

Next, we compute the cross-correlation amplitudes of 
the VERITAS data with these AT SZ and T syn maps. On 
average, the VERITAS data show no correlation with the 
synchrotron map, and an insignificant, — 1.0a (—1.6a) 
anti-correlation with the SZ map (i.e. a positive corre- 
lation with the y-paramctcr). Here and below, results 
in parenthesis are the conservative confidence levels after 
correcting for the wobble-mode background subtraction, 
as described in §4. 

Features seen in both SZ and synchrotron maps are 
more pronounced in the Western part of the cluster, pos- 
sibly due to some Eastern foreground (e.g., see Figure 
5a) or CMB fluctuation. Consequently, for RA < 195.3 
(no more than 0.5° East of the cluster's center), the SZ 
anti-correlation strengthens to —1.6a (—2.3a). Note that 
this signal does not arise from the spatial separation be- 
tween the central SZ decrement and the peripheral 7-ray 
ring, as this would correspond to a positive correlation 
here. 

More importantly, when we split the data into two re- 
gions, outside and inside the inner b — 1° ellipse bound- 



ing the 7-ray signal, correlations emerge between 7-rays 
and both radio signals. Outside this inner ellipse, i.e. 
along the 7-ray ring, the VERITAS map correlates with 
the synchrotron map at the +2.8a (+3.9ct) confidence 
level, and anti-correlates with the SZ map at the — 2.6a 
(—3.6a) level. The signal is dominated by the Western 
half of the ring, where these correlations strengthen to 
+3.2<7 (+4.1a) with the synchrotron and — 3.9ct (—5. Oct) 
with the SZ. The results depend somewhat on the precise 
region examined. For example, including the 0.5° region 
east of the center, the synchrotron correlation strength- 
ens to +3.5ct (+4.6ct) while the SZ correlation remains 
unchanged. 

Within the b = 1° ellipse, the VERITAS map docs 
not correlate with the SZ signal, but does show a 
— 1.9a (—2.7a) anti-correlation with the synchrotron 
map, mainly in the Western part, where it reaches —2.2a 
(—2.8a). The inner anti-correlation between the VER- 
ITAS and synchrotron maps, evident in the 7-ray un- 
dcrluminous radio halo and radio relic, and their outer 
positive correlation, arc illustrated in Figure 9. 

In the periphery of the cluster, the positive correla- 
tion between the VERITAS map and the synchrotron 
signal agrees with predictions for the coincident inverse- 
Compton and synchrotron emission from the same rela- 
tivistic, shock-accelerated electrons, as we show in §8 and 
§9 below. The anti-correlation with the SZ map agrees 
with the anticipated SZ signal from the thermal gas, ex- 
pected to cutoff sharply as the thermal pressure drops 
beyond the virial shock, as shown in §9. 

In the inner parts of the cluster, an anti-correlation be- 
tween 7-ray and synchrotron radiation was anticipated 
in secondary electron models, where electrons are con- 
tinuously injected into the ICM in both strongly and 
weakly magnetized regions 22 . Here, electrons magne- 
tized above (below) £> cm b = (Snucmb) 1 ^ 2 ~ 3 /xG, where 
« cm b is the CMB energy density, lose most of their en- 
ergy to synchrotron (inverse-Compton) radiation in the 
radio (7-ray) band 26 , leading to alternating radio- and 
7-ray-bright regions, as discussed in §10. 



8. THE OBSERVED 7-RAY SIGNAL AGREES 
WITH PREDICTIONS FOR AN £ e m ~ 5% 
ELECTRON INJECTION RATE 



Cluster virial shocks last for long, ~ H^ 1 Hubble 
timcscales, and are strong, in particular where cold gas 
accretes from the voids and the Mach number squared is 
3> 10. Therefore, such shocks are thought to accelerate a 
flat, dn/de oc e -2 spectrum (constant energy per logarith- 
mic energy interval) of rclativistic electrons up to high, 
cooling-limited energies. The dominant cooling process 
here is Compton scattering off CMB photons, implying 
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FIG. 8: Radio maps of the Coma region. 

the low- frequency AT SZ = — 2y in mK) and synchrotron 
T syll (i//23 GHz) 3 2 in mK) maps, spectrally extracted from the WMAP seven-year all sky map- L . 
showing the 7-ray bright region, are the same as in Figure la, and the cyan contours show the VERITAS mosaic region. 



Thermal SZ (left; colorbar is 



(right; colorbar is 
The dashed curves, 



an electron energy cutoff 1 

^ T + 1 ^ / 3eBk B T 

2 y (T - l)mpa T u cmb 

~ 60(T 10 S .i) 1/2 TeV , (5) 

where k B T = 10T 10 kcV and B = O.LB .i fiG arc the 
temperature and magnetic field amplitude downstream 
of the shock. Here, e and m e are the electron charge 
and mass, gt is the Thompson cross section, and c is 
the sound velocity. The photon energy E ~ 220 GeV 
which maximizes the VERITAS sensitivity corresponds 1 1 
to electrons with energy e ~ (i?/3fcBT cm b) 1 ^ 2 TO e c 2 ~ 
9 TeV, below but not far from the estimated cutoff. 

The main uncertainty in modelling the nonthermal 
emission from a virial shock stems from the presently 
poor understanding of particle acceleration and magneti- 
zation in collisionless shocks. The inverse-Compton emis- 
sion from the shock can be estimated by assuming that 
a fraction £ e = 0.05£ e ,5 of the post-shock thermal energy 
is deposited in relativistic electrons. Electrons with en- 
ergies ^ 100 MeV cool quickly, in a narrow shell behind 
the shock, so the logarithmic emissivity per unit surface 
area of the shock is 

E 2 d 3 N _ 3!j e vnk B T 
dtdAdE ~ 4\n(e m ax/m e c 2 ) 

~ lO-^sn-slfo 72 erg s" 1 cm" 2 , 

where n = 10~ 5 n_5 cm -3 is the (thermal) electron num- 
ber density, and v = (fcsT/12/i) 1 / 2 is the downstream 



velocity. 

The 7-ray signal from a Coma-like cluster can be com- 
puted using Eq. (6) in the framework of a cosmological 
simulation ' , without introducing additional free param- 
eters. However, the comparison between simulated and 
observed clusters is complicated by a secondary uncer- 
tainty, involving the local accretion rate of the cluster. 
This rate varies spatially along the shock front, fluctu- 
ates in time, and differs even among similar clusters as a 
function of their specific environment. We parameterize 
the average accretion rate using the dimensionless vari- 
able 

. _ M 47rR 2 ^nv 
m = ~MH ~ f b MH ' (7) 

where here we assumed a spherical shock of radius R. 
In the j3- model discussed in Appendix A, m = 1.4, and 
in the simulated cluster' rh ~ 0.3. These are averaged 
quantities, with strong local fluctuations seen in simula- 
tions and expected in practice. 

Consider a planar cut parallel to the line of sight, along 
which the virial shock is approximately a ring of radius R 
and of negligible thickness. A beam that intersects the 
cut at distances {ri,?^} < R from the center covers a 
shock segment of length 2R\ arcsin(r 2 /i?) — arcsin(ri //?) | . 
Therefore, at the near edge of a long cylindrical shock of 
angular radius ip, a beam of area 8 2 would detect > E 
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FIG. 9: 

Composite 7-ray and synchrotron emission. 

Superimposed VERITAS map (green; colorbar range: 
— 3.4ct to +3. Oct) and synchrotron map (red; spectrally 
extracted from the WMAP 7-ycar full sky-map; col- 
orbar range: —1.3 to +5.4(^/GHz) 32 K). showing a 
positive correlation (yellow) outside the inner edge of 
the 7-ray ring (dashed contour) and an anti-correlation 
(spatially separated green and red, the latter domi- 
nated by the radio halo and relic) inside the ring. 



energy photons with number flux 

J(> E ) c 7r - 2aTCSin{1 - 9M E 



d 3 N 
dt dAdE 



(8) 



4tt(1 + z) 4 9/ip 

A similar result is expected in the prolate shock inferred 
in Coma, where for the VERITAS map 

/ / \ 1/2 

-8 ( 



J(> 220 GcV) ~ 2 x 10- 



\ioe 



xTfg 72 s" 1 cm -2 ster -1 



(9) 

In the isothermal /3-model, this may be written as J(> 
220 GeV) ~ 5 x 10" 8 £ e , 5 m s" 1 cm" 2 ster" 1 . 

In order to compare this estimate with the 7-ray struc- 
ture seen in the VERITAS map, recall that a Ict fluctua- 
tion corresponds to an excess of N a = — 32 photons 
with E > 220 GeV in a beam of radius 6 = 0.2°. The 
map was produced by a t = 18.6 hour exposure 1 , with 
an A ~ 4 x 10 8 cm 2 effective area (at ~ 220 GeV, see 39 ). 
A Ict signal in the beam therefore corresponds to a pho- 
ton flux 

J la ~ ~ 3.1 x 10" 8 s" 1 cm" 2 ster" 1 . (10) 

TTV^tA 

The average signal inferred from the histogram (Figure 
4), fJ-2 ~ Ml — 1.4ct, is thus comparable to the predicted 
signal in Eq. (9) for standard parameters. 



The agreement between the predicted (Eq. 9) and ob- 
served (Eq. 10) 7-ray flux, and the similar significance of 
the 7-ray features in the observed (Figure la) and simu- 
lated (Figure lb) maps, suggest that to within an uncer- 
tainty factor of a few, £ e m ~ 5%. Localized regions with 
~ 3ct detection, and possibly brighter emission masked by 
the wobble-mode background subtraction, suggest that 
at least locally the product ^ e m can exceed its inferred 
average value by a factor of a few. 

The total flux from the virial shock can be estimated 
by taking 9 — > 2ip in Eq. (8), and integrating over the 
solid angle enclosed by the shock. For a cylindrical ring 
of radius b and length £ b, this gives 



F(>E) 



2b 2 C d 3 N 
:E- 



4(1 + z) 4 dtdAdE 

2Q-9 C^l^e^-^T^ o _ 1 ^-2 



(11) 



E^l + zY 

where bx = 6/1° and E x = E/l GcV. For the VERITAS 
band and b = 1.3°, we use the isothermal /3-model to ob- 
tain F(> 220 GeV) ~ 3 x 10" n C£ e , 5 77i s" 1 cm" 2 . Cali- 
brating the free parameters against the local Ict VER- 
ITAS flux (Eq. 10), this becomes F{> 220 GeV) ~ 
2 x 10" n C s" 1 cm" 2 . 

Extrapolating the VERITAS-calibrated signal to the 
EGRET band (assuming a flat spectrum), yields F(> 
100 MeV) ~ 4 x 10" 8 C s" 1 cm" 2 , comparable to the 
4 x 10" 8 s" 1 cm" 2 EGRET upper limit on emission 
from Coma 10 ' 11 . The corresponding F(> 200 MeV) ~ 
2 x 10 _8 C s" 1 cm" 2 exceeds the 5 x 10~ 9 s" 1 cm" 2 up- 
per limit imposed by Fermi 12 , but this limit was obtained 
using an aperture optimized for emission from the cen- 
tral, X-ray bright region of the cluster. Note that at such 
low energies, the Galactic foreground is not negligible, 
and both instruments lack the high angular resolution 
of VERITAS essential in order to identify the ring and 
avoid foreground contamination. 

The above analysis assumed a flat spectrum of accel- 
erated electrons, leading to a flat (p = 2) photon spec- 
trum. Indeed, comparing the VERITAS measurement 
with the EGRET upper limit implies that the photon 
spectral index between 100 MeV and 220 GcV cannot 
be much softer than p = 2.1. As the VERITAS mea- 
surement probes only photons near 220 GeV, the above 
estimates strictly pertain to electrons in energies near 
10 TeV. Here, we find a logarithmic electron injection 
efficiency 



: de 



ln(e maa; / trnin) 



0.3% 



(12) 



9. THE SYNCHROTRON AND SZ SIGNALS 
AGREE WITH PREDICTIONS FOR £s ~ 1% 
MAGNETIZATION 

The same electrons that scatter CMB photons up to 
the 7-ray band also emit synchrotron radiation, as they 
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gyrate in the shock amplified magnetic fields. The ra- 
tio between radio and 7-ray emission from a strong virial 
shock is approximately given by the ratio between mag- 
netic and CMB energy densities, 



vL, 



B 2 /8ir 
EJ(> E) ~ 

^cmb 



(13) 



In order to compute the synchrotron signal, some as- 
sumption must be made concerning the magnetic field 
behind the shock. Typically 4,1 , one assumes that a frac- 
tion £b = 0.01£_b,i of the thermal energy is deposited 
in downstream magnetic fields, based on observations of 
supernova remnant shocks; some magnetization is also 
needed to accelerate the 7-ray emitting electrons. How- 
ever, this assumption was not tested so far near the virial 
shock, and the value of £b is not well constrained. 

We may use the coincident 7-ray and radio measure- 
ments to estimate the magnetic field at the shock. The 
strongest signal is found in the lowest WMAP frequency, 
v = 23 GHz, where bright spots along the shock show 
T syn ~ 0.1 mK, corresponding to a brightness 



ul v = 2— ksTgyn ~ 10 9 erg s 1 cm 2 ster" 



(14) 



The same bright spots are detected by VERITAS 
at ~ 3a, corresponding to EJ(> E) ~ 3.2 x 
10 -8 erg s -1 cm -2 ster -1 . Equation (13) then yields 
B ~ 0.4 /xG, which according to the /3-model corresponds 
to £ B ~ 0.5%. 

We may use the SZ signal near the shock to test the 
above estimates. For example, at v = 23 GHz one ex- 
pects the ratio between the synchrotron emission and 
the SZ decrement to be T syn /AT SZ ~ — 0.3^,5^8,1^ for 
an isothermal sphere distribution 1 '', and T syn /AT SZ ~ 
— 0.5£e :5 £B,i for a /3-model (see Appendix A). 

Inspecting the coincident synchrotron (T syn ~ 0.1 mK 
at 23 GHz) and SZ (AT SZ ~ —0.2 mK) signals measured 
along the 7-ray ring, we find that AT syn /AT sz ~ —0.5. 
This supports the above estimates £ e m ~ 5%, and 
£b — !%• Note that the parameter £b refers to the 
magnetic energy weighted over the relativistic electron 
distribution, and so is elevated by clumping effects : . 



10. IN THE INNER REGIONS, THE 
7-RAY— SYNCHROTRON ANTI-CORRELATION 
SUPPORTS STRONG MAGNETIZATION AND 
HADRONIC INJECTION 

Within 6=1°, we find an anti-correlation between 7- 
rays and synchrotron emission. This signal is not highly 
significant (—2.2a and — 2.8a, before and after correcting 
for wobble background removal), so deeper observations 
are needed to clarify its nature. Here we consider the 
implications of the anti-correlation assuming it is indeed 
real. In such a case, as the synchrotron signal involved 
is dominated by the radio halo and radio relic, some of 



the VERITAS signal must arise from the inner ICM, well 
inside the virial radius. 

In the radio band. Coma is a classical example of a 
merging cluster, harboring both an extended radio halo 
and a peripheral radio relic 1 ' 1 . The > Mpc size of such ra- 
dio halos implies an ongoing injection of relativistic elec- 
trons, but it is debated whether they arise from hadronic 
collisions involving CR ions (secondary models 22 ' 25 ' 26 ' 40 ) 
or from in-situ (re)accclcration in weak shocks or tur- 
bulence (primary models 21 ' 24 ' 41 ). Secondary models can 
also explain radio relics - " and minihalos"' as arising from 
the same distribution of CR protons, if the latter are 
homogeneously distributed 22 . Here, radio-bright regions 
correspond to strong, B > £? cm b magnetic fields. The ex- 
istence of such strongly magnetized regions is also being 
debated, in Coma 12 ' 45 and in other clusters 22 ' 44-46 . 

In secondary models, electron injection is an ongoing 
process, proportional to the gas density but largely in- 
dependent of shocks, turbulence, and other changes in 
the flow. Therefore, the sum of synchrotron and inverse- 
Compton emission from the ICM varies slowly across the 
cluster, whereas their ratio can show strong variations, 
proportional to (B / B cm ^) 2 , leading to a radio-7-ray anti- 
correlation. In contrast, in primary models, the electrons 
are injected only locally, in shocks or turbulence where 
the magnetic field is amplified, and cool long before the 
magnetic field can decay. Hence, a positive correlation 
is expected here between synchrotron and 7-ray signals; 
inverse- Compton emission far from the radio sources is 
expected to be negligible. 

In both secondary and primary models, 7-ray emission 
is expected also from 7r° decay, following CR ion colli- 
sions with the ambient nucleons. This signal is expected 
to be smoother than the inverse-Compton emission, and 
so more susceptible to wobble-mode background sub- 
traction, in particular if the CRs are distributed ho- 
mogeneously, as inferred from hadronic radio models 22 . 
For example, in a /3-model with homogeneous CR ions, 
wobble-mode removes > 50% of the signal in the central 
0.2° beam, which corresponds to the magnetized core, 
and > 90% of the signal beyond tp = 0.3°. The tt° sig- 
nal, independent of the local magnetic field and stronger 
in the center, should positively correlate with the syn- 
chrotron map, which is dominated by the central radio 
halo. It cannot account for the anti-correlation presently 
observed. 

Therefore, the inner 7-ray— synchrotron anti- 
correlation, if confirmed, supports the notion of 
magnetized, > 3 /iG regions in the ICM, and favors 
a secondary electron injection model. Here, the pres- 
ence of a weak 7-ray signal far from the center (but 
inward of the virial shock) supports a homogeneous CR 
proton distribution. In non-magnetized regions, the 
T syn (23 GHz) ~ 0.2 mK synchrotron signals from the 
radio halo and relic translate to an inverse-Compton flux 
J(> 220 GeV) ~ 3x 10 -9 s -1 cm -2 ster -1 . By Eq. (10), 
this corresponds to an Sj ~ O.ltr VERITAS significance 
per 8 = 0.2° beam, so 7-ray emission around the core 
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can plausibly lead to the observed anti-correlation. 

Integrated within a b = 1°, £ = 3 ellipse, emission from 
secondary electrons corresponds to an S ~ 0.4ct nominal 
VERITAS signal, assuming homogeneous CR ions and a 
/3-model gas distribution. Such a signal can be further 
tested in a deeper 7-ray observation. If wobble-mode 
corrections can be minimized, the (~ 4 times stronger ' ) 
7T° signal may be measured as well. 



11. SUMMARY AND DISCUSSION 

We analyze the VERITAS mosaic (Figure la) of the 
Coma cluster, and find a large-scale, extended 7-ray fea- 
ture around the cluster. It is best described as an elon- 
gated (C > 3), thick (Ab > 0.5°) elliptical ring, with 
semi-minor axis b ~ 1.3° (Figure 2), oriented approxi- 
mately along the East- West direction (4> ~ —5°; Figure 
3), towards the LSS filament connecting Coma with Abell 
1367 (Figure 6a). The 7-ray ring is seen at a nominal 
3.3<7 confidence level, but comparison with mock VERI- 
TAS noise maps implies a higher, 4.5ct confidence. Fur- 
ther correcting for excessive wobble-mode background 
subtraction, we conclude that the ring presents with a 
confidence level of at least 5.7a. If we do not assume 
a-priori that the ring is aligned with the LSS filament, 
the significance becomes 4.2ct (> 5.4a corrected). Our 
wobble-mode corrections, derived in §3 and §4 (see Fig- 
ure 4, and Eq. 3), are quite conservative. 

The size and morphology of the 7-ray ring agree with 
predictions for the signature of the virial shock, produced 
as shock-accelerated electrons inverse Compton-scatter 
CMB photons. The VERITAS map thus correlates with 
a normalized simulated 7-ray map of a Coma-like cluster 
(Figure lb), once the main LSS filaments are similarly 
oriented (Figure 6), at a 3.7a (4.9a) confidence (no free 
parameters; see Figure 7). The best fit is obtained for 
ellipticity £ ~ 4.5, but we cannot rule out a cylindrical 
shock parallel to the LSS filament. The brightness of the 
signal agrees with predictions as well, provided that the 
rate of relativistic electron energy deposition is £ e m ~ 5% 
of the thermal energy per Hubble time (within a factor 
of a few) . 

Interestingly, a 3a signal found 1 ' by stacking the 
EGRET data around 447 rich clusters corresponds, 
when compared with source number counts computed 
analytically 1 and using cosmological simulations 1 , to an 
average £ e ,5W ~ 4%, in agreement with the value we infer 
in Coma. 

Inverse-Compton emission from virial shocks can ex- 
plain the hard X-ray signals observed in several clusters 1 ', 
provided that £ e < 10%. For example, £ e m ~ 4% 
was recently inferred from hard X-ray emission from 
a LSS filament 20 . The VERITAS ring corresponds to 
a ~ 5C cV s" 1 cm -2 flux at 20-80 keV energies (as- 
suming a flat spectrum), accounting for the (8.1 ± 
2.5) eV s _1 cm -2 hard X-ray signal observed 19 . High 
resolution observations in hard X-rays, in particular with 



NuSTAR, should resolve the virial ring in Coma, and de- 
tect virial shocks in many other clusters. 

An analysis of the WMAP seven-year data reveals 
a positive correlation between the 7-ray ring and syn- 
chrotron emission (Figures 8b and 9), at the 2.8ct (3.9ct) 
confidence level, and an anti-correlation with the SZ map 
(Figure 8a) at — 2.6er (— 3.6ct). The radio signals are 
clearer in the Western half of the ring, towards the LSS 
filament, where the correlations reach 3.5ct (4.6a) and 
—3.9a (—5. Oct), respectively. They agree with predic- 
tions for the synchrotron emission from the 7-ray emit- 
ting electrons and for the SZ decline with thermal pres- 
sure beyond the shock, if the downstream magnetic field 
is B ~ 0.4 fiG, corresponding to £b ~ 0.5% magnetiza- 
tion (within a factor of a few). 

The SZ signal we find, extending in some parts out to 
ip ~ 2.3° away from the center, agrees with the SZ detec- 
tions by WMAP 21 and Planck 1 '. However, the periph- 
eral synchrotron signal coincident with the 7-rays and 
with the SZ cutoff was not reported so far. The 7-ray 
ring and these radio features are much farther out than 
the ift ~ 0.5° weak shocks seen in X-ray, radio 1 ' 1 , and 
SZ 17 data. 

Inward of the ring, the VERITAS map shows a 
— 1.9ct (— 2.7a) anti-correlation with the synchrotron 
map, stronger in the Western region where C = —2.2a 
(—2.8a). This is evident in the radio halo and radio relic, 
which appear undcrluminous in 7-rays (Figure 9). Such a 
signal was predicted in models where relativistic electrons 
are injected into the ICM continuously, in both strongly 
and weakly magnetized regions which consequently ap- 
pear radio bright/7-ray faint and vice versa. While the 
present signal is weak and requires confirmation, it is con- 
sistent with analytic estimates. The extent of the anti- 
correlation and its magnitude support secondary electron 
injection by collisions of cosmic-ray protons with the am- 
bient nucleons, and favor a homogeneous CR distribu- 
tion. 

A VERITAS observation deeper than the present 18.6 
hour exposure would decisively test our results, in partic- 
ular if wobble-mode corrections can be avoided or min- 
imized, for example by using a larger, w > 1° wobble 
radius. This would allow for a more precise reconstruc- 
tion of the virial shock structure and a determination of 
variations in the electron deposition rate along the shock 
front. Such data would be valuable in the study of large- 
scale structure formation at low redshift, mapping the 
feeding pattern of Coma, probing the surrounding voids, 
the Coma-Abell 1367 filament, and other LSS filaments, 
and tracing the warm-hot intcrgalactic medium (WHIM) 
immediately behind the shock. 

Such a deep observation would test the anti-correlation 
we find between 7-ray and radio maps at small, b < 1° 
radii. A significant anti-correlation would confirm the 
hadronic, secondary electron model for cluster radio 
sources, and provide a direct measure of magnetic varia- 
tions in the ICM. A deep observation would thus probe 
the CR proton distribution out to large radii, address the 
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origin of these protons, and in particular place an upper 
limit on proton acceleration in the virial shock. 

Observations at lower energies, for example using the 
~ 50 GeV MAGIC Cerenkov telescopes 1 *, should be 
able to utilize the higher photon flux to clearly identify 
the 7-ray ring. At yet lower energies, the Fermi 7-ray 
space telescope is marginally sensitive to the signal at 
< 1 GeV, but it may be challenging to resolve the shock 
structure with the > 1° (68% containment) point spread 
function 49 . 

The VERITAS energy threshold, E ~ 220 GeV, is 
not far below the expected photon cutoff, E max ~ 
lOTio-Bo.i TeV. Observations at higher energies would 
become increasingly sensitive to this cutoff, thus probing 
the conditions at the shock and the nature of particle 
acceleration, although pair production off the infrared 
background also becomes important at such energies'". 
Interestingly, no evidence for extended emission was 
found in a 1.1 TeV observation of Coma by the HESS 
telescopes 28 . Despite the higher energy threshold and 
the shorter, 8.2 hour observation, it was only ~ 2.3 times 
less sensitive than the VERITAS observation (assuming 
a flat spectrum), thanks to the larger effective area 51 , 
~2x 10 9 cm 2 . The HESS map shows no significant sig- 
nal or correlation with the VERITAS data, but the 7-ray 
ring may well be shallowly buried under the noise. Note 
that the optical depth of the HESS photons due to pair 
production off the infrared background, estimated 50 as 
t ~ 0.2, may not be negligible. 



Finally, consider the contribution of cluster virial 
shocks to the cxtragalactic 7-ray and radio backgrounds. 
These components were computed analytically 1-3 , 
and calibrated numerically using cosmological 
simulations ' ' . While the main model parameters £ e m 
and £b are likely to fluctuate within clusters and vary 
among different clusters, it is instructive to adopt the val- 
ues inferred in Coma as typical. An accretion rate m ~ 1 
is intermediate between the analytic model 1 and the low 
redshift result of numerical simulations . Normalizing 
all clusters by this value, we obtain a diffuse 7-ray com- 
ponent E 2 dJ/dE ~ 0.3(£ e m/5%) keV s" 1 cm" 2 ster" 1 , 
dominating the extragalactic background . In 
the radio, the revised parametrization yields a 
vJ„~3x 10- 11 (^to/5%)(^ b /1%) erg s" 1 cm" 2 stcr" 1 
synchrotron signal, dominating the extragalactic 
low frequency radio background ' . It is observable 
through 6T1 ~ l(^/GHz)~ 3 K fluctuations at multipoles 
400 < I < 2000 with present interferometers such as 
LOFAR and EVLA. 
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Appendix A: /? model for the shock 



An analytic model for emission from the virial shock requires some assumptions specifying the gas distribution. 
Simple choices include an isothermal sphere 1 or a isothermal /3- model 47 , the latter of which we discuss below. 
In our simple model, it is assumed that (i)the electron number density is approximately given by 



r 2 



n(r) = n ^1 + (Al) 

out to the virial radius; and (ii) on average, the gas distribution is approximately isothermal and in hydrostatic 
equilibrium. In Coma 52 , the central density is approximately no = (3.0 ± 0.6) x 10~ 3 cm -3 , the core radius is 
r c = 343^20 kpc, the index /3 = 0.654to'o2i' an d the (emission measure- weighted) temperature is T = 8.38±0.34 keV. 

With the above assumptions, the region enclosing k = 200^200 times the critical mass density p c is characterized 
by 

' ~ 2.3«at» 2 M P C . ( A2 ) 



;/3Tr 



M K ~ ~ 1.4 x 1O 15 K2O 1 O /2 M , (A3) 



and 

n K = n(r K ) ~ 7.2 x 10~ 5 k 2 oo cm -3 . (A4) 

We may now compute the accretion parameter, 

_ 4irrlpn K v K 1/2 
mK ~ f b M K H - L4K 200: (A5) 
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the inverse- Compton emissivity per unit shock area, 



E 2 d 3 N _ 3M b mM K Hk B T 
dtdAdE I6irr 2 /.iln(e max /m e c 2 ) 



4.6 x 10 8 ?ii£ ei 5K2Qo er g s 1 cm 2 



and the magnetic field at the shock, 



B K = \IZ B 8^nk B T ~ 0.6^ 2 ^i M« . (A7) 
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